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Abstract 

MicroCraft  Inc.  performed  a  specialized  calibra¬ 
tion  of  the  NASA  Langley  104B  balance  (fabricated 
from  CVM  200  alloy)  In  their  Automatic  Balance 
Calibration  System  (ABCS)  machine  to  identify  and 
determine  the  magnitude  of  anelastic  effects  in  the 
calibration  results.  The  results  analyzed  thus  far 
are  for  the  axial-force  gage.  These  results  show 
that  anelastic  effects  for  the  balance  in  question 
are  substantial  enough  to  consider  revising  calibra¬ 
tion  methodology  for  situations  where  highest 
accuracy  is  needed.  Worst  case  effects  found  thus 
far  are  of  the  order  of  0.1  percent  balance  gage 
capacity.  Modeling  of  the  anelastic  relaxation  pro¬ 
cess  with  time  constants  that  are  linear  with  stress 
level  was  successful  in  reducing  the  uncertainty  by 
at  least  50-percent. 


Introduction 


cent  of  capacity.  Although  it  may  actually  be  worse 
because  of  bias  that  results  from  noncalibrated  fac¬ 
tors  such  as  temperature,  anelastic  effects  (hyster¬ 
esis),  and  inertial  effects.  Opportunity  for  improve¬ 
ment  includes  calibration  for  these  effects,  new 
materials,  and  new  methods  of  sensing  forces  and 
moments  such  as  optical-based  methods  (birefrin- 
gent  materials  and  doped  optical  fibers  with  either 
Bragg,  or  long-wave  grating.)  This  paper  reports  on 
progress  in  accounting  for  anelastic  effects  in  the 
relationship  between  applied  forces  and  moments 
and  strain-based  voltage  output.  Steinle^  examined 
the  application  of  the  fundamental  relationship 
involving  stress,  strain,  and  their  respective  time- 
derivatives  presented  by  Zener^  to  the  anelastic 
calibration  of  the  axial-force  gage  of  a  1.1-in.,  13-8 
PH  stainless  steel  balance.  Zener’s  relation  is: 
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The  wind  tunnel  testing  community  has  gener¬ 
ally  accepted  the  desirability  of  being  able  to  mea¬ 
sure  forces  and  moments  accurately  enough  to 
support  definition  of  drag  coefficient  accurate  to 
within  one  count  (0.0001).  The  two  dominant  items 
affecting  the  ability  to  measure  drag  with  sufficient 
accuracy  is  the  measurement  of  angle  of  attack  to 
better  than  0.01  deg  and  axial  force  to  better  than 
0.1  percent  of  rated  capacity  of  the  balance.  Gains 
in  accuracy  of  both  measurements  are  being 
sought.  Optical  methods  for  angle-of-attack  mea¬ 
surement  are  coming  into  use,  as  are  corrections 
to  accelerometer-based  angle  measurements  for 
inertial  effects.  As  for  balances,  uncertainty  is  gen¬ 
erally  thought  to  be  of  the  order  of  0.2  to  0.25  per- 


Here,  o,  e,  and  C  represent  stress,  strain,  and 
modulus,  respectively,  and  and  K2  are  time  con¬ 
stants.  Zener  used  “C”  for  modulus  Instead  of  the 
traditional  “E”  for  Young’s  modulus.  Zener’s  con¬ 
vention  is  retained  since  it  may  prove  necessary  to 
model  “C”  as  a  function  of  strain  and/or  tempera¬ 
ture  in  the  future  to  reach  a  more  robust  model  of 
the  anelastic  process.  Equation  (1)  with  “C”  being  a 
constant  can  be  solved  in  integral  form  for  either 
stress,  given  strain  as  a  function  of  time,  or  vice- 
versa.  This  was  done  in  Ref.  1,  with  stress 
assumed  proportional  to  force  and  strain  propor¬ 
tional  to  output  voltage.  Equation  (2)  is  the  result¬ 
ing  integral  expression  for  strain  as  a  function  of 
stress: 


*  The  research  reported  herein  was  performed  by  the  Arnold  Engineering  Development  Center  (AEDC),  Air  Force  Materiel 
Command.  Work  and  analysis  for  this  research  were  performed  by  personnel  of  Sverdrup  Technology,  Inc.,  AEDC  Group,  technical 
services  contractor  for  AEDC.  Further  reproduction  is  authorized  to  satisfy  needs  of  the  U.  S.  Government. 
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If  the  time  base  (tq)  is  started  from  the  point 
where  stress  can  be  approximated  by  a  constant 
value,  then  Equation  (3)  simplifies  to: 

(t-tg) 
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where  with  initiai  values  of  e  and  a  =  0,  can  be 
found  from  Equation  1 .1  by  integrating  from  0  to  tq 
to  result  in: 
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From  Eqs.  (2)  and  (3),  it  is  seen  that  with 
Zener’s  formulation,  the  amount  of  anelastic  strain 
can  be  normalized  by  the  amount  of  stress.  This 
means  that,  regardless  of  the  amount  of  stress, 
strain-based  instrumentation  will  thus  experience 
the  same  percentage  of  error,  independent  of  the 
working  stress  level  of  the  balance.  Material  prop¬ 
erties  that  result  in  the  time  constants  are  the  key 
to  reduction  of  anelastic  effects. 


This  present  study  investigated  the  magnitude 
of  anelastic  effects  that  might  be  present  in  results 
from  a  typical  longitudinal  aerodynamic  test  series 
performed  in  the  AEDC  16T  wind  tunnel  and 
applied  the  method  to  reduce  anelastic  error.  The 
balance  used  is  the  NASA  Langley  104B  balance. 
Figure  1  is  a  sketch  of  the  104B  Balance.  Details  of 
the  balance.  The  approach  taken  in  this  study  was 
to  perform  a  self-consistent  calibration  of  the  104B 
balance  in  the  MicroCraft  Force  Measurement  Sys¬ 
tems  ABCS  machine.  The  study  was  accomplished 
by  first  checking  the  machine’s  load  cells  for 
anelastic  effects  and  then  performing  a  series  of 
calibrations  of  the  balance.  This  series  started  with 
a  standard  calibration,  followed  by  a  time-corre¬ 
lated  calibration  of  the  balance  that  followed  a  load 
profile  over  the  time  period  representative  of  an 
entire  start-run-stop  series  from  a  test  of  a  low 
aspect-ratio  aircraft  configuration  at  high  subsonic 
speeds.  Unlike  the  wind  tunnel  load  versus  time 
profile,  this  series  of  loads  was  applied  at  approxi¬ 
mately  constant  intervals  of  time.  Three  time  peri¬ 
ods  were  employed.  Two  of  these  time  periods 
produced  a  total  time  interval  representative  of  the 
wind  tunnel  test  (90  and  99  minutes.)  The  third 
time  period  (40  minutes)  resulted  from  the  ABCS 
operating  at  its  standard  rate.  The  series  finished 
with  a  check  of  the  relaxation  characteristics  of  the 
104B  balance  at  constant  load  with  both  load 
applied  and  with  load  removed. 


In  Ref.  1 ,  the  relaxation  process  taking  place  in 
the  1.1 -in.  balance  for  constant  load  was  best 
described  by  two  modes  (two  time  constants 
required  instead  of  just  K2)  that  differed  by  roughly 
an  order  of  magnitude  (nominally  100  sec  and 
1000  sec.)  Zener’s  equation  was  easily  modified  to 
produce  an  additional  relaxation  term.  The  reason 
for  the  two-mode  representation  was  not  identifi¬ 
able.  Time  constants  for  the  magnitude  of  the  load 
(absolute  value)  increasing  were  found  to  be 
somewhat  different  for  magnitude  reducing.  The 
results  showed  that  the  hysteresis  loop  uncertainty 
could  be  reduced  by  over  84  percent.  The  reader  is 
referred  to  Ref.  1  for  details.  These  results  were 
quite  encouraging.  However,  further  development 
was  clearly  indicated  before  a  viable  calibration 
technique  for  a  six-component  balance  could 
emerge  from  this  line  of  work. 


NASA  LaRC  Balance  NTF-104B 

The  NASA  LaRC  NTF-104B  was  constructed 
for  use  in  Langley’s  National  Transonic  Facility  in 
1 993.  The  design  is  based  on  LaRC  standard  prac¬ 
tices  outlined  in  Ref.  3.  An  outline  drawing  that  pro¬ 
vides  the  balance’s  physical  characteristics,  minus 
strain  gages  and  wiring,  is  shown  in  Fig.  1 .  In  addi¬ 
tion,  the  balance  full-scale  load  capacities,  com¬ 
bined,  are  listed  in  Fig.  1 . 

The  balance  material  is  VascoMax  200  type  C 
that  was  heat-treated  to  a  Rockwell  C  of  42-45. 
This  material  is  utilized  for  all  of  the  NTF  series  bal¬ 
ances  because  of  its  toughness  characteristics  at 
cryogenic  temperatures.  The  balance-to-model 
interface  is  a  cylindrical  fit  that  must  duplicate  the 
fit  of  a  ring  and  plug  gage  set.  A  dowel  pin,  typically 
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General  Note: 

All  dimensions  are  in  inches 


Figure  1.  NASA  Langley  104B  balance. 

a  0.0002-in.  interference  fit,  secures  the  balance  to 
the  model  and  defines  its  location.  The  sting  end 
attachment  is  tapered  with  a  key  and  set  screw 
flats  for  positioning  and  securing. 

SK  350-ohm  type  strain  gages  from  Micro  Mea¬ 
surements  Group,  along  with  their  M610  adhesive, 
are  utilized  on  the  balance.  Current  procedures  for 
installing  strain  gages  on  cryogenic  balances  are 
detailed  in  Ref.  4  (balance  NTF-104B  varies  from 
these  procedures  in  some  areas  due  to  the  time  it 
was  strain  gaged  in  1993).  One  variation  is  the 
moisture-protective  coating.  The  104B  has  a 
Teflon-based  coating  that  has  been  repaired  with 
Mcoat  B,  a  Micro  Measurements  Group  Nitril  rub¬ 
ber  based  compound.  (Note:  NTF-104B  is  in  need 
of  being  refurbished  to  NASA  LaRC’s  new  strain 
gaging  standards  developed  since  1993.)  The  bal¬ 
ance  also  contains  three  temperature  sensors, 
thermocouples.  They  are  located  at  three  stations 
along  the  balance  length,  forward,  middle,  and 
rear.  These  are  used  to  monitor  balance  tempera- 


Reference  5  discusses  LaRC’s  calibration 
methods. 

Anelastic  Calibration  Test  Description 
MicroCraft  ABCS  Machine  &  Load  Ceiis 

220 

The  ABCS  is  the  primary  calibration  system 
for  internal  balances  at  MicroCraft  Force  Mea- 
surement  Systems  in  San  Diego,  CA.  The  ABCS 
(Fig.  2)  is  a  fully  self  contained  and  integrated  6- 
deg-of-freedom  calibration  system.  Figure  2a  is  a 
photograph  of  the  overall  ABCS,  and  Fig.  2b  is  a 
close-up  view  of  the  load  adapter  rig.  The  load 
cells  used  to  sense  applied  loads  are  clearly  visible. 
The  system  is  capable  of  performing  all  tasks 
required  for  the  balance  calibration  process.  The 
machine  applies  a  force  to  the  balance  using 
hydraulic  actuators.  Precision  load  cells  and  high- 
resolution  optical  sensors  are  used  to  determine 
applied  load  vector  magnitude  and  direction.  Load 


a.  Overall  system 


tures  and  gradients  during  calibration  and  wind 
tunnel  testing. 

Calibration  of  the  NTF-104B  in  LaRC’s  manual 
calibration  stand  in  1993  generated  the  following 
accuracy  results  at  room-temperature  conditions 
(two-standard  deviation  accuracy  was  determined 
from  back  calculated  loads  versus  applied  loads). 


Normal 

0.14% 

Axial 

0.23% 

Pitch 

0.15% 

Roll 

0.16% 

Yaw 

0.11% 

Side 

0.28% 

b.  Load  adapter 

Figure  2.  Automatic  Balance  Calibration  System. 
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position  is  determined  from  mechanical  measure¬ 
ments.  The  ABCS  is  based  on  a  non-repositioning 
principle.  The  applied  force  actuators  are  allowed  to 
move  as  a  result  of  balance  deflection.  The  position 
of  the  deflected  balance  is  measured  relative  to  the 
fixed  end  of  the  balance  axis  using  six,  0.5-micron 
incremental  optical  linear  displacement  gages.  By 
measuring  this  deflection,  the  applied  force  vectors 
may  be  established  mathematically.  The  ABCS  is 
essential  to  the  full  characterization  and  calibration 
of  multicomponent  balances.  The  ABCS  can  apply 
single  and  multicomponent  loads  in  any  combina¬ 
tion  within  a  balance’s  defined  load  limits. 

Load  Cell  Calibration  Procedure 

Prior  to  using  the  load  cells  in  the  ABCS  sys¬ 
tem,  each  cell  was  calibrated  to  determine  primary 
sensitivities,  linearity,  and  anelastic  characteristics. 
Since  the  anelastic  characteristics  are  time  based, 
each  load  cell  was  loaded  to  its  maximum  as 
quickly  as  possible.  Once  full-scale  load  was 
reached,  load  cell  output  was  monitored  at  speci¬ 
fied  time  intervals.  Table  I  shows  the  data  interval 
schedule.  In  all,  data  were  recorded  over  a  30-min 
window  at  full  load.  Once  completed,  the  weight 
was  removed  from  each  load  cell  and  the  same 
data  interval  schedule  was  repeated  at  zero  load 
conditions  (recording  relaxation).  This  process  was 
repeated  for  both  tension  and  compression  on 
each  load  cell. 


Table  I  -  Data  Interval  for  Load  Cell  Calibration 


Beginning  Time 

Ending  Time 

Data  Interval 

0 

1  minute 

5  seconds 

1  minute 

5  minutes 

30  seconds 

5  minutes 

10  minutes 

1  minute 

10  minutes 

30  minutes 

5  minute 

Balance  Calibration  Process 

The  calibration  process  was  accomplished 
using  both  the  ABCS  and  manual  calibration  pro¬ 
cedures.  Balance  matrix  coefficients  were  deter¬ 
mined  for  a  variety  of  matrix  models  using  both 
ABCS  and  manual  data.  The  anelastic  static  load¬ 
ing  was  completed  using  dead  weight  for  the  appli¬ 
cation  of  a  three-  component  load  over  the  same 
time  variable  used  during  the  load  cell  calibrations. 


The  ABCS  calibration  process  included  combi¬ 
nations  of  two-component  loadings  to  full-scale 
balance  loads.  Both  a  6x27  and  a  6x96  calibration 
matrix  were  derived  from  this  calibration  data.  A 
second  calibration,  also  using  two  component 
loads,  was  completed  to  a  reduced  loading  range 
designed  to  match  the  maximum  loads  typically 
achieved  during  wind  tunnel  testing.  The  balance 
was  also  loaded  to  full  scale  using  an  837-point 
loading  schedule.  This  load  schedule  was  selected 
as  representative  of  the  standard  calibration  that  is 
preformed  on  this  balance.  The  837-point  schedule 
applies  loads  in  a  different  combination  to  define 
the  balance’s  characteristics.  The  balance  was 
then  loaded  using  a  calibration  schedule  represen¬ 
tative  of  actual  load  conditions  recorded  during  a 
wind  tunnel  test.  The  loads  were  applied  over  the 
same  time  as  that  of  the  test.  Five  complete  cali¬ 
brations  of  the  test  condition-loading  schedules 
were  performed,  with  the  fourth  schedule  run  at  the 
standard  speed  setting  on  the  ABCS.  The  balance 
was  also  loaded  manually  with  dead  weights  to 
verify  the  matrix  produced  by  the  ABCS  calibration. 
Over  the  course  of  three  days,  a  total  of  6200  load¬ 
ing  points  were  applied  to  the  baiance. 

Analysis  of  Test  Results 

Load  Cell  Anelastic-Relaxation  Time-Constant 
Calibration 


If  a  linear  relation  between  stress  in  the  load 
cell  and  applied  load  is  assumed  along  with  a  lin¬ 
ear  relation  between  strain  sensed  by  the  bridge 
system  of  the  load  cell  and  output  voltage.  Equa¬ 
tion  (3)  for  constant  load  can  be  rewritten  as: 


V  =  KF. 


(x-tn) 


1-e 


+  V,  e 


Ko 


(5) 


where  V  represents  the  normalized  output  voltage 
(millivolts/volt),  K  is  the  cell  calibration  constant 
relating  force  to  millivolts/volt,  and  F  represents  the 
applied  calibration  force.  The  constants  K,  K2,  and 
Wxq  were  determined  from  a  least-squares  fit  of  the 
test  results.  The  change  in  output  due  to  anelastic 
effects  is  then  KF^^j  -  ,  Application  of  Eq.  (5) 

required  interpretation  to  remove  outliers  attributed 
to  noise  in  the  output.  Weight  pan  motion  and 
ground  vibration  and/or  air  pressure  variation 


4 

American  Institute  of  Aeronautics  and  Astronautics 


AIAA-99-0682 


resulting  from  commercial  aircraft  operations 
nearby  inducing  motion  in  the  building’s  calibration 
support  structure  are  thought  to  have  been  a  factor 
in  the  noise  that  seemed  to  be  present  in  the  test 
resuits.  The  caiibration  test  was  not  equipped  to 
detect  such  accelerations  by  independent  means.  It 
should  be  noted  that  a  milli-g  of  vertical  acceleration 
would  induce  a  noise  spike  in  the  test  results  of 
approximately  0.1  percent  of  applied  ioad.  This 
value  of  noise  is  an  order  of  magnitude  above  the 
floor  noise  level  sought.  With  these  caveats,  it 
appeared  that  the  load  cells  used  in  the  ABS 
machine  are  very  good  and  have  hysteresis  that  is 
not  statisticaiiy  significant.  Results  from  the  load 
cell  calibration  were  therefore  not  used  in  the  anai- 
ysis  and  are  not  given  herein. 

Balance  Anelastic  Relaxation  Time-Constant 
Calibration 

Difficulties  with  the  method  used  to  determine 
the  reiaxation  time  constants  are  associated  with 
noise  in  the  data  (predominantiy  inertial  effects), 
change  in  applied  force  due  to  relaxation,  and 
interaction  factors  due  to  multicomponent  balance 
loading.  The  procedures  employed  to  deal  with 
these  difficulties  required  developing  a  model  of 
the  process  and  then  utilizing  regression  to  deter¬ 
mine  constants  in  the  model.  The  methodology 
used  to  compensate  for  the  above  difficulties  is  as 
follows: 

Inertial  Effects.  As  above,  judgment  was 
employed  to  remove  what  was  considered  spuri¬ 
ous  data  by  treating  them  as  outliers  and  removing 
them  through  regression  anaiysis.  To  accomplish 
removai  of  outiiers,  a  ieast-squares  fit  of  the  error 
process  was  used  to  provide  an  initial  estimate  that 
was  then  manually  examined  to  identify  outliers, 
and  a  final  least-squares  fit  of  the  results  was 
obtained. 

Change  In  Applied  Force.  Relaxation  effects 
changing  the  applied  load  vector  during  the  anelas¬ 
tic  relaxation  manual  loading  are  modeied. 
Because  of  the  ioading  approach  taken  which  pro¬ 
duced  a  combination  of  axiai-force,  normal-force, 
and  pitching-moment  ioading,  relaxation  in  pitch 
angle  results  in  a  time  variation  of  appiied  load.  A 
variation  of  angle  within  0.01  deg  pitch  angle 


change  could  not  be  sensed  with  the  available 
instrumentation.  Since  defiection  change  due  to 
reiaxation  was  expected  to  be  less  than  0.01  deg, 
output  of  the  angle-of-attack  sensor  was  not  moni¬ 
tored.  Determination  of  pitch  angie  to  the  order  of 
0.001  deg  is  needed.  The  total  elastic  defiection  of 
the  balance  due  to  the  loading  is  estimated  at  0.4 
deg  (based  on  deflection  information  from  the 
ABCS  caiibration.)  Assuming  an  anelastic  effect  of 
the  order  0.1  percent  of  balance  capacity,  this 
would  lead  to  a  change  in  angle  of  0.0012  deg, 
which  is  not  detectable  through  either  the  angle 
sensor  used  or  regression  analysis  and  hence,  is 
not  considered  important  in  the  analysis  of  results. 
Nevertheless,  it  was  modeled,  checked  by  regres¬ 
sion,  and  found  to  not  be  statistically  significant. 

Time  Constant  Determination.  Assuming  that 
relaxation  in  pitch  angle  is  negligible,  only  relax¬ 
ation  at  constant  stress  in  the  balance  flexures  is 
required.  Since  the  balance  is  homogenous,  the 
time  constant  associated  with  all  gages  of  the  bal¬ 
ance  is  the  same  for  the  same  conditions  of  strain 
and  temperature.  This  also  includes  pitch  angle 
relaxation  since  the  balance  elastic  properties  that 
produce  strain  also  result  in  rotation.  The  relax¬ 
ation  is  expected  to  be  different  for  load  being 
applied  as  opposed  to  load  being  removed.  Exami¬ 
nation  of  the  anelastic  loading  results  indicates  that 
this  is  the  case.  Inertial  effects  appear  to  dominate 
the  axial-force  relaxation  at  constant  load;  hence,  it 
was  not  possible  to  extract  a  meaningful  time  con¬ 
stant  from  the  axial-force  load  relaxation  with  load 
applied.  However,  the  normal-force  relaxation 
results  were  affected  much  less,  and  use  of  Eq.  (5) 
resulted  in  a  time  constant  value  of  nominally  565 
sec  for  applied  load.  A  time  constant  of  480  sec 
produces  a  reasonable  fit  of  the  data  for  both  the 
normal-force  and  axial-force  gages  for  relaxation 
after  load  is  removed.  The  relaxation  fit  of  the  data 
for  both  gages  is  shown  in  Fig.  3.  All  data  points 
taken  are  included.  The  difference  in  character 
between  the  axial-force  loading  shown  in  Fig.  3a, 
as  compared  to  the  normal-force  loading  in  Fig.  3b 
can  be  seen.  The  time  constant  numbers  indicated 
should  not  be  construed  as  absolute  since  the 
noise  in  the  data  makes  it  impossible  to  determine 
the  time  constant  within  a  close  tolerance.  A  better 
set  of  laboratory  conditions  and  additional  instru¬ 
mentation  is  expected  to  produce  better  results. 
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Figure  3.  relaxation  time  constant  determination. 

However,  the  results  indicated  are  quite  sufficient 
for  the  purpose  intended. 

There  is  an  inherent  assumption  in  this  process 
that  the  relaxation  constants  are  the  same,  regard¬ 
less  of  the  strain  (or  stress)  level.  If  the  constants 
are  not,  then  that  variation  must  be  modeled  suffi¬ 
ciently  well,  and  the  above  results  only  give  a  value 
appropriate  to  the  strain  state  of  the  material.  Fur¬ 
thermore,  the  relaxation  constants  probably  are  a 


function  of  temperature  as  well.  The  experiment 
reported  herein  was  at  constant  temperature  and 
thus  could  not  provide  any  information  suggestive 
of  the  significance  of  temperature  on  relaxation.  As 
will  be  seen  later,  modeling  of  the  relaxation  pro¬ 
cess  under  general  loading  has  led  the  authors  to 
the  conclusion  that  a  higher-order  modeling  of 
relaxation  constant  behavior  as  a  function  of  stress 
is  needed  to  match  the  results  of  this  investigation. 
A  linear  model  of  time  constant  versus  strain  (or 
stress)  was  proven  to  be  necessary,  but  not  suffi¬ 
cient  to  capture  the  entire  nature  of  observed 
results. 


Investigation  of  Criticality  of  Load  vs.  Time. 

Application  of  the  same  load  schedule  utilizing  suf¬ 
ficiently  different  time  intervals  indicative  of  the  dif¬ 
ference  between  calibration  machine  and  wind  tun¬ 
nel  data  acquisition  rates  will  show  whether  or  not 
time  is  a  significant  factor  in  balance  calibration. 
The  procedure  followed  was  to  perform  the  “same” 
ABCS  calibration  load  schedule  with  two  different 
time  intervals  between  each  new  load.  The  com¬ 
plete  load  series  is  representative  of  the  entire  load 
cycle  experienced  in  a  wind  tunnel  test  cycle  from 
start  of  the  tunnel  to  stop  with  the  104B  balance. 
Figure  4  shows  the  representative  load  conditions 
for  normal  force,  pitching  moment,  and  axial  force. 
The  load  variations  are  typical  of  a  transonic  air¬ 
craft  test  that  encompasses  early  compressible 
flow  with  strong  leading-edge  suction  (Polar  E) 
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Figure  4.  Nominal  load  schedule. 
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through  buffet  conditions  (Polar  I).  The  range  of 
pitch  plane  loading  conditions  approximated  a  third 
of  the  capacity  of  the  balance.  The  remaining  three 
load  components  were  small  and  are  not  shown. 
Unavoidably,  there  were  small  differences  between 
the  six  components  of  load  being  applied  at  com¬ 
parable  conditions.  Thus,  it  was  necessary  to  cor¬ 
rect  the  output  voltage  for  the  difference  in  load 
schedule.  This  was  done  by  multiplying  the  inverse 
of  a  6X6  calibration  matrix  times  the  difference  in 
loading  to  arrive  at  a  voltage  increment  to  apply  to 
the  output  generated  at  the  highest  rate.  It  is 
assumed  that  after  this  correction  procedure,  tak¬ 
ing  the  difference  between  the  two  output  voltages 
will  remove  bias  in  the  calibration  results,  leaving 
the  relative  anelastic  effects. 


A  comparison  of  the  error  in  calibrated  axial 
force  (applied-computed  load)  is  shown  in  Fig.  5 
for  the  ABCS  standard  rate  (40  min  total  for  408 
data  points)  against  the  maximum  W.T.  load  cycle 
(99  min,  total).  Output  and  load  cycle  output  (40 
min,  total),  corrected  to  the  W.T.  load  cycle,  is 
shown  in  Fig.  5  for  three  of  the  nine  polars.  Absent 
anelastic  effects,  the  dominant  error  should  be  bias 
in  the  calibration  that  would  result  in  all  of  the  data 
scattered  about  a  line  of  agreement  with  45-deg 
slope.  As  can  be  seen,  this  does  not  occur.  There 
is  considerable  difference  for  each  of  the  polars, 
which  is  attributed  to  anelastic  effects. 
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Figure  5.  Effect  of  different  load  application  rates 
on  axial  force  calibration  error. 


The  error  result  for  each  load  rate  contains  the 
random  error  difference,  plus  calibration  bias  error, 
plus  the  difference  in  anelastic  effects  between  the 
two  load  versus  time  histories.  The  range  of  rela¬ 
tive  anelastic  error  is  substantial.  The  peak  differ¬ 
ence  is  about  0.4  lb  which  is  0.13  percent  of  gage 
capacity.  Since  the  peak  loading  was  about  28-per- 
cent  capacity,  0.4  lb  difference  is  thus  equivalent  to 
0.46  percent  of  peak  load.  Consequently,  the  vari¬ 


ability  between  the  two  load  rates  is  0.23  percent  of 
peak  load,  which  is  significant,  even  though  it  is 
only  0.065  percent  of  gage  capacity.  Hence,  any 
notion  of  bias  error  in  the  balance  calibration 
results  that  is  based  on  calibration  performed  at  a 
high  productivity  that  does  not  match  the  wind  tun¬ 
nel  load  profile  productivity  could  be  substantially 
unrealistic,  unless  the  anelastic  effect  can  be  mod¬ 
eled  and  both  the  calibration  and  test  results  can 
be  corrected.  If  modeling  is  not  practical,  then  the 
obvious  approach  where  higher-than-customary 
accuracy  is  needed  would  be  to  actually  “fly”  the 
balance  in  an  automatic  calibration  machine  that 
has  the  capability  to  vary  load  versus  time  in  a 
manner  that  will  closely  follow  the  balance  output 
voltage  time-history.  The  current  software  control¬ 
ling  the  ABCS  machine  does  not  permit  this  type  of 
profile.  Conceivably,  it  is  possible  to  do  this  with 
current  modern  control  technology  that  has  learn¬ 
ing  ability. 

Modeling  Anelasticity  in  Calibration  Resuits. 
Assuming  that  Eq.  (1)  essentially  represents  the 
physics  involved  in  modeling  anelasticity,  then  it  is 
obvious  that  some  sort  of  continuous  numerical  or 
analog  integration  process  will  be  required  to  cor¬ 
rect  wind  tunnel  data  for  anelasticity.  In  either 
case,  the  constants  (or  variables)  and  K2  must 
be  determined. 


Determination  of  Stress-Rate  Constant  (K.|). 
In  addition  to  the  relaxation  time-constant  for  con¬ 
stant  load  (K2),  it  is  necessary  to  determine  the 
stress  rate  constant  to  be  able  to  represent  the 
anelastic  effects.  Two  different  rates  of  application 
of  load  are  required  to  determine  this  time  con¬ 
stant.  By  modeling  the  load  versus  time  process, 
the  calibration  results  discussed  above  were  used 
regressively  to  determine  the  stress-rate  constant. 
The  approach  taken  was  to  assume  that  the  stress 
at  any  location  in  the  balance  Oj  is  linearly  related 
to  the  six  components  of  forces  and  moments,  Fj 
as 


6  3ai 
=  1  ' 


(6) 


Then  this  result  is  substituted  into  a  more  con¬ 
venient  form  of  Eq.  (2),  along  with  the  assumption 
that  output  voltage  is  proportional  to  strain.  Eq.  (2) 
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is  integrated  by  parts  to  produce  a  result  conve¬ 
nient  to  discretizing  the  load-time  history: 

f  (^-y^ 


T  Xf\ 


e 


V 


Letting  superscript  t  represent  time  and  sub¬ 
script  j  the  voltage  from  gage  j  at  time  x,  and  sub¬ 
stituting  Eq.  (6)  into  Eq.  (7)  along  with  the  assump¬ 
tion  that  each  contribution  to  stress  produces  a 
corresponding  contribution  to  voltage  through  the 
induced  strain,  Eq.  (8)  is  obtained: 


6  6 


(^-y 


(t-y 


avp 


Ko 


+  V 

j 


K, 


'-iqO' 

,  i  =  li  V  iA  2J 


(8) 


The  partial  derivatives  of  output  voltage  with 
respect  to  force  or  moment  in  Eq.  (8)  are  approxi¬ 
mated  by  the  linear  terms  of  the  inverse  calibration 
matrix.  Should  there  be  significant  cross-product 
terms,  these  terms  would  also  be  added.  Eq.  (8)  is 
in  a  form  that  can  be  converted  to  finite  differences. 
Assuming  that  over  each  interval  of  time  the  rate  of 
change  of  each  force  can  be  treated  as  a  constant, 
Eq.  (8)  can  be  integrated: 
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Equation  (9)  shows  that  the  anelastic  effect  is 
embodied  in  the  second  through  the  last  term  on  the 
right-hand  side.  Taking  the  limit  as  T/K2  becomes 
very  large  at  constant  F,  Eq.  (9)  reduces  to  just  the 


first  term,  which  is  the  true  “back-calibration”  result 
with  anelastic  effects  removed.  Consequently,  dur¬ 
ing  the  calibration  process,  a  subtraction  of  the  volt¬ 
age  bias  due  to  anelasticity  would  produce  a  data 
set  that  can  be  used  to  arrive  at  the  correct  value  for 
the  balance  calibration  matrix  with  anelastic  effects 
removed.  Owing  to  the  formulation,  the  predicted 
anelastic  component  is  coupled  to  the  current  esti¬ 
mate  of  the  end  result.  Consequently,  rigorous  pro¬ 
cessing  of  the  calibration  results  requires  iterating 
the  recomputed  calibration  results  back  into  Eq.  (9). 
Iteration  should  not  be  necessary,  since  the  magni¬ 
tude  of  the  anelastic  component  is  of  the  order  of 
0.1  percent  of  gage  capacity.  Use  of  Eq.  (9) 
requires  testing  to  determine  whether  or  not  the 
magnitude  of  the  voltage  is  greater  or  less  than  the 
fully-relaxed  value.  If  it  is  greater,  then  the  K2  value 
corresponding  to  decreasing  load  (e.g.,  480  sec  at 
zero  load)  is  required.  If  it  is  less,  then  K2  equals  the 
value  corresponding  to  increased  loading. 


Since  the  true  value  of  voltage  with  anelastic 
effects  cannot  be  known  unless  is  also  known, 
subtracting  the  result  of  two  identical  load  sched¬ 
ules  performed  at  different  time  intervals  will  thus 
eliminate  the  first  term  on  the  right-hand  side  of  Eq. 
(9),  including  bias.  A  least-squares  approach  can 
then  be  taken  to  determine  K^,  which  will  best 
match  the  differential  anelastic  voltage  inferred 
from  Fig.  5.  The  approach  taken  is  to  start  from  the 
condition  of  all  forces  and  voltage  being  zero  and 
then  marching  forward  in  time  for  each  data  point. 


The  assumption  of  constant  values  for  and 
K2  did  not  prove  satisfactory  herein.  Assuming  that 
the  general  form  of  Eq.  (1)  is  a  good  representa¬ 
tion,  the  next  logical  approach  is  to  treat  and  K2 
as  variables.  A  linear  variation  with  absolute  value 
of  strain  (or  stress)  for  K-j  and  K2  was  assumed: 


(10) 


^2”  ^02^"^  ^02^^  +^2^1) 

The  use  of  the  absolute  value  was  chosen 
because  the  output  of  each  balance  gage  is  from  a 
full  bridge  that  is  subjected  to  tensile  stress  in  half 
of  the  legs  of  the  bridge  and  compressive  stress  in 
the  other  half.  Thus,  the  change  in  the  output 
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should  not  vary  with  the  sign  of  the  stress  level.  It 
should  vary  with  whether  or  not  relaxation  is 
increasing  or  decreasing  the  absolute  value  of 
stress.  The  rate  of  change  of  the  relaxation  con¬ 
stants  is  slow  enough  that  they  can  be  treated  as 
constants  over  each  interval  of  time  and  updated 
accordingly.  This  was  accomplished  by  substituting 
Eqs.  (10a)  and  (103b)  into  Eq.  (9).  Regression  to 
determine  the  best  overall  fit  of  all  of  the  anelastic 
computed  results  against  the  difference  in  output 
for  the  two  experiment  load  rates  was  performed. 
The  values  arrived  at  are  shown  in  Table  II. 


Table  II.  Anelastic  Time  Constants 


Constant 

Koi 

Ci 

Ko2 

^2 

Increasing  Load 

565.9 

0.424 

564.1 

0.424 

Decreasing  Load 

480 

0.395 

480.1 

0.4398 

As  previously  indicated,  the  corrected  result 
was  used  in  the  regression  to  determine  the  above 
final  set  of  constants.  Figure  6  shows  the  compari¬ 
son  of  the  experimentally  determined  difference  in 
anelastic  results,  along  with  the  computational 
results  based  on  the  above  constants  for  the  axial- 
force  gage  for  all  nine  polars.  The  polars  (A 
through  I)  are  in  the  order  in  which  the  load  was 
applied.  In  general,  it  is  obvious  that  the  approach 
taken  is  in  the  correct  direction,  and  it  does  show 
promise  in  that  the  effect  is  predicted  to  a  large 
degree.  There  are  instances,  e.g.,  polars  A,  D,  and 
E,  where  the  trends  predicted  are  significantly  dif¬ 
ferent  from  the  experimental  results.  This  argues 
for  more  definitive  work,  especially  in  determining 
material  properties  from  direct  measurement. 

Figure  7  shows  the  relative  anelastic  error  for 
both  experiment  and  computed  results  sequentially 
for  all  polars.  The  agreement  with  the  general  trend 
shows  that  in  the  gross  sense,  the  model  captures 
the  trend.  A  fine  detail  comparison  between  experi¬ 
ment  and  computed  results  shows  that  a  different 
model  is  required  to  properly  capture  the  true 
nature  of  the  anelastic  effects. 

A  projection  of  the  absolute  level  of  anelasticity 
for  each  of  the  two  loading  rates  is  shown  in  Fig.  8 
for  polars  E,  F,  and  I.  The  projection  is  based  on 
the  constants  given  in  Table  II.  The  difference 
between  the  two  sets  of  data  is  then  the  computed 
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Figure  6.  Comparison  of  computed  on  experimen¬ 
tal  anelastic  error  due  to  load  rate. 
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difference  in  aneiasticity  for  the  two  ioad  rates.  The 
maximum  value  of  predicted  anelastic  error  is  0.01 
millivolts/volt,  which  is  equivalent  to  0.19  percent 
gage  capacity!  Further  work  is  needed  to  extract  a 
meaningful  estimate  of  the  absolute  level  of 
aneiasticity  in  each  of  the  load  conditions.  One 
impediment  is  the  lack  of  definitive  knowledge  of 
the  timing  for  each  data  point.  The  analysis  was 
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Figure  6.  Concluded. 


performed  assuming  that  each  point  was  taken  at 
the  same  interval.  A  further  difficulty  is  associated 
with  the  need  to  do  a  more  careful  analysis  of  the 
calibration  results  and  extract  a  first  estimate  of  the 
anelastic  effects  from  the  input  data  and  then 
reprocess  the  calibration  to  try  to  remove  bias  error 
that  is  not  due  to  aneiasticity.  This  will  be  tedious  at 
best.  On  the  other  hand,  the  usefulness  of  the 
present  method  for  correcting  for  aneiasticity  can 


Figure  7.  Sequential  relative  anelastic  error. 
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Figure  8.  Predicted  anelastic  error. 
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be  seen  by  comparing  Figs.  9  and  10.  The  reduc¬ 
tion  in  relative  anelastic  error  between  the  two  load 
application  rates  is  seen  to  be  at  least  50  percent. 
Although  not  conclusive,  it  does  offer  promise  of 
obtaining  similar  results  on  an  absolute  basis. 
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Figure  9.  Experimental  relative  anelastic  error. 
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Figure  10.  Residual  experimental  relative  anelas¬ 
tic  error  after  computational  correction. 

Concluding  Remarks  and  Recommendations 

1.  It  is  hoped  that  the  result  of  this  work  will 
spark  interest  in  others  to  join  in  researching  the 
correction  for  anelasticity,  which  shows  promise  of 
reducing  the  uncertainty  associated  with  strain- 
based  measurements  of  forces  and  moments.  As  a 
first  step,  it  is  recommended  that  a  relaxation  check 
be  made  for  balance  output  under  load  and  with 
load  removed  to  determine  both  the  magnitude  of 
relaxation  and  the  time  constants  involved  and  that 
any  results  of  such  tests  be  shared  with  the  AIAA 
GTTC,  Balance  Calibration  Working  Group. 


shown  herein.  Considerably  more  work  is  war¬ 
ranted,  especially  direct  laboratory  measurement 
of  anelastic  characteristics  under  the  range  of  con¬ 
ditions  experienced  by  a  balance. 

3.  The  mean  value  of  anelasticity  at  each  data 
point  in  the  calibration  introduces  a  bias  at  each 
load  in  the  calibration  process  that  then  results  in 
a  nonlinear  term  being  in  the  gage  calibration 
when  there  would  otherwise  not  be  one.  Removal 
of  anelastic  effects  will  then  make  it  possible  to 
establish  the  true  nonlinearity  in  the  calibration, 
which  could  result  from  a  number  of  causes, 
including  the  modulus  not  being  a  constant  and 
the  inherent  nonlinearity  in  a  basic  strain  gauge 
output.  In  either  case,  the  current  calibration  meth¬ 
odology  employed  to  process  the  post-anelastic 
correction  calibration  results  will  handle  those  sys¬ 
tematic  nonlinear  effects. 

4.  Although  this  paper  has  not  discussed  appli¬ 
cation  of  the  anelastic  calibration  results  to  the 
reduction  of  wind  tunnel  data,  it  does  not  pose  a 
problem.  The  basic  differential  equation  still  holds, 
and  it  is  trivial  to  process  the  results  going  from 
measured  voltage  to  a  corrected  voltage  that  can 
then  be  used  with  a  standard  calibration  to  arrive 
at  the  final  result. 

5.  The  success  of  the  present  method  in  reduc¬ 
ing  apparent  anelastic  error  between  two  different 
calibration  load  time  cycles  merits  further  work  to 
develop  a  calibration  technique  and  the  neces¬ 
sary  application  of  anelastic  correction  to  test 
results  and  to  transition  such  a  technique  to  rou¬ 
tine  application.  Any  efforts  toward  this  effect 
should  be  shared  with  the  AIAA  GTTC  Balance 
Calibration  Working  Group. 
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